† Corresponding author. E-mail:
Project supported by the National Natural Science Foundation of China (Grant No. 11574083).
In this work, we theoretically probe into the photo-induced hydrogen bonding effects between S0 state and S1 state as well as the excited state intramolecular proton transfer (ESIPT) behavior for a novel 2-[1,3]dithian-2-yl-6-(7 aH-indol-2-yl)-phenol (DIP) probe system. We first study the ground-state hydrogen bonding O–H⋯N behavior for DIP. Then we analyze the primary geometrical parameters (i.e., bond length, bond angle, and infrared (IR) stretching vibrational mode) involved in hydrogen bond, and confirm that the O–H⋯N of DIP should be strengthened in the first excited state. It is the significant prerequisite for ESIPT reaction. Combining the frontier molecular orbitals (MOs) with vertical excitation analyses, the intramolecular charge transfer (ICT) phenomenon can be found for the DIP system, which reveals that the charge redistribution facilitates ESIPT behavior. By constructing potential energy curves for DIP along the ESIPT reactional orientation, we obtain quite a small energy barrier (3.33 kcal/mol) and affirmed that the DIP molecule undergoes ultrafast ESIPT process once it is excited to the S1 state and quickly transfers its proton, forming DIP-keto tautomer. That is why no fluorescence of DIP can be observed in experiment, which further reveals the ultrafast ESIPT mechanism proposed in this work.
Hydrogen bonds are ubiquitous and essential in the field of chemistry and biochemistry, and govern chemical reactions, supramolecular molecules, molecular assemblies, and life process.[1–3] The strength of hydrogen bonding depends on the electronegativity of donor atom which the hydrogen atom is connected with, and also depends on the ability to absorb the atom which the hydrogen atom interacts with through hydrogen bond. Both intramolecular and intermolecular hydrogen bonding affect the absorption and emission properties of the molecules. Often the fluorescent based probe molecule which forms intramolecular and intermolecular hydrogen bond is sensitive to the external environment. Therefore, there has become an increasing interest in generating novel hydrogen bonding sites in molecular probe and assemblies.[4–8]
As one of the most elementary processes throughout nature, proton transfer (PT) belongs to a fundamental class of chemical field taking place along with intramolecular or intermolecular hydrogen bond.[9–14] In view of the photo-excitation behavior, the photo-induced PT (i.e., excited state intramolecular proton transfer (ESIPT)) has been regarded as one of most significant photochemical and photo-physical and photo-biological reactions in a variety of reactions.[15–20] Therefore, it has become one of the hotspots of experimental and theoretical researches in the field of chemistry. Due to the ground-state transient property, molecules containing ESIPT characteristics have been adopted in a variety of applications. The basic photo-physical and photo-chemical role of ESIPT is that the proton of one existing hydrogen bond transfers from the donor atom to acceptor moiety, which facilitates forming a transitory equilibration between enol configuration and keto configuration, emitting double fluorescence peaks and large Stokes shift.[21–30]
Generally, enol configuration should be photo-excited to S1-state enol* structure. Subsequently, the ESIPT happens and forms keto* structure in the S1 state. In view of the differences between enol* configuration and keto* configuration, normally, keto* configuration emits a longer wavelength fluorescence.[31–38] Therefore, the dual emission band could be detected. It covers a broad region of steady-state emission by promoting ESIPT molecules to become suitable for cell images, UV filters, white emitting OLEDs, fluorescent probes, molecular switchers, etc.[21–38] Therefore, it cannot be denied that dynamics of ESIPT reaction is beneficial for more novel applications in future.
Taking advantage of the modulations of fluorescence as well as the sensitivity to surrounding environment for ESIPT materials, advanced applications have been made in various novel areas as mentioned above. Recently, Chang et al. reported a new mercury ion probe based on the ESIPT mode.[39] Based on the classical 2-(2-hydroxyphenyl)benzoxazole and 2-(2-hydroxyphenyl)benzothiazole molecules and through the two-step reaction, the novel 2-[1,3]dithian-2-yl-6-(7 aH-indol-2-yl)-phenol (DIP) probe can be synthesized.[39] They inferred that the probe response might be due to the de-protection of dithian function in aqueous solution. Also, using smartphone as a signal-capturing device, the detection of metal ions is tested based on the DIP system. Thus, it can be seen that DIP molecule owns good prospects and potential for more novel development of probes. Before developing new purposes, the most fundamental excited state dynamical behaviors of DIP itself should be the primary one to be explored. In other words, a clear dynamical process of excited state for the DIP system could push and facilitate novel applications. To the best of our knowledge, only spectroscopic techniques containing time-resolved as well as steady-state electronic spectra could just reveal oblique aspects about molecular photochemical properties.[21–38] Up to now, previous investigations and observed phenomenon of DIP system are limited except experiment study. In view of the significance of theoretical work, we adopt the quantum chemical computational manner to study photo-induced cases relevant to hydrogen bonding interactions as well as ESIPT reaction. The ground-state investigations are based on density functional theory (DFT), and excited-state explorations are implemented by time-dependent DFT (TDDFT) method. In this work, we mainly focus our attention on configurations and chemical bond changes involved in intramolecular hydrogen bonding interactions. Moreover, all the chemical vibrations and infrared (IR) vibrational spectra, vertical excitation energy, frontier molecular orbitals (MOs), charge density and electronic density redistribution, and related ESIPT kinetic curves of DIP system are calculated to analyze the straightforward information about excited state dynamical process and ESIPT mechanism. Our simulated results show an ultrafast ESIPT behavior for DIP, which explains experimental fluorescence phenomenon.
In the present work, all the calculations about electronic structure are carried out by using the Gaussian 09 program suit.[40] The geometric optimization of DIP-enol and DIP-keto systemsare performed in the S0 state by using DFT and in the S1 state via TDDFT method. Here, we should mention that the reason why we select the TDDFT method is because it has become a very useful method to investigate the hydrogen-bonding interactions occurring in the excited states. Becke’s three-parameter hybrid exchange functional with the Lee–Yang–Parr gradient-corrected correlation (B3LYP functional) and the triple-ζ valence quality with one set of polarization function (TZVP) basis set are selected in both DFT process and TDDFT process.[41–43] Any constrain is imposed on none of symmetry, bond lengths, bond angles or dihedral angles in the geometric optimization calculations. In order to evaluate the solvent effects and make them consistent with previous experimental results, acetonitrile solvent is selected through our calculations based on the polarizable continuum model (PCM) by using the integral equation formalism variant (IEFPCM).[44–46] It is worth mentioning that the IEFPCM model mitigates the computational burden of explicitly modeling solvent molecules and the specific interactions with the solute. The solvent can be treated as a constant dielectric reaction field where the charge density of the solute is projected onto a grid on the surface of a solvent cavity and polarized, depending on the value of the solvent dielectric. The resulting polarized charges on the cavity affect the charge density of the molecule, etc., until self-consistency is achieved.
All the local minima are confirmed by the absence of an imaginary mode in the vibrational analysis calculations. In addition, the S0- and S1-state potential energy curves are scanned by constrained optimization and frequency analyses to obtain the thermodynamic corrections in their corresponding electronic states by keeping the O–H bond length in step of 0.05 Å. Zero-point energy corrections and thermal corrections to the Gibbs free energy are also carried out according to the harmonic vibrational frequencies. Harmonic vibrational frequencies in the S0 and S1 states are determined via the diagonalization of Hessian matrix. The excited-stated Hessian matrix is obtained by the numerical differentiation of analytical gradient via central differences and default displacement of 0.02 Bohr. The infrared intensity is determined from the gradient of the dipole moment. Fine quadrature grids of size 4 are employed. The self-consistent field (SCF) convergence thresholds of the energy for both the ground state and excited state optimization are set to be 10−8 (default settings are 10−6). Harmonic vibrational frequencies in the ground and excited state are determined by diagonalizing the Hessian matrix.
As shown in Fig.
In order to better explore the excited state dynamical behavior and ESIPT process, we first pay attention to the excited state hydrogen bonding behaviors for the DIP-enol molecule. As mentioned in previous classical work, the obtained hydrogen bond lengths, bond angles and other geometrical parameters,[48–57] and the corresponding optimized chemical bond parameters involved in intramolecular hydrogen bond O–H⋯N of DIP-enol compound are all listed in Table
![]() | Table 1. Simulated bond lengths and bond angles of DIP-enol and DIP-keto structure in S0 and S1 states, respectively. . |
Monitoring the variations about vibrational spectral shifts is also an effective method to study the hydrogen bonding dynamics. Particularly, the vibrational spectral shifts involved in specific vibrational mode can also reflect the chemical bonding changes in excited states.[48–57] Therefore, the infrared (IR) spectra of DIP-enol compound in vibrational region of O–H stretching mode are also calculated and displayed in Fig.
![]() | Fig. 4. Theoretical IR vibrational spectra for (a) O–H vibrational mode of DIP-enol structure and for (b) H–N vibrational mode of DIP-keto structure, in both S0 state and S1 state. |
When we mention the photo-excitation process, we also theoretically predict the low level absorbing transitions via TDDFT/B3LYP/TZVP (i.e, vertical excitation energy of S0–S1, S0–S2, …, S0–S9) in acetonitrile solvent. Given the deactivations in high excitations, we only focus on the first excited state transition. Listed in Table
![]() | Table 2. Calculated electronic excitation energy (in unit nm), corresponding oscillator strengths and the relative compositions of the low-lying three excited states for compound DIP. . |
In order to clearly understand the ESIPT mechanism, the potential energy kinetic curves with O–H⋯N are scanned via optimization manner. To be specific, the scanning of ESIPT curves is through constrained optimization by keeping O–H bond length. Our scanned potential energy curves are shown in Fig.
In this work, we carefully explore photo-induced hydrogen bonding dynamics and ESIPT behavior for a DIP system. The dynamics about excited state hydrogen bonding effects mainly focuses on the geometries and IR vibrational spectra. Changes of chemical parameters reveal that the hydrogen bond O–H⋯N of DIP should be strengthened via photo-excitation, which should be the precondition of ESIPT reaction. Furthermore, the photo-induced charge redistribution and ICT properties promote the ESIPT tendency. Particularly, the increased electronic densities around proton acceptor N moiety facilitate the absorption of hydrogen proton. The analyses about potential energy curves along the possible ESIPT path demonstrate that the ESIPT behavior of DIP system should be ultrafast due to low barrier. According to the ultrafast ESIPT mechanism, we can reasonably explain previous experiments. This work clarifies the excited state dynamical behavior for DIP system, but also fills the vacancy and inadequacy of previous experiment.
[1] | |
[2] | |
[3] | |
[4] | |
[5] | |
[6] | |
[7] | |
[8] | |
[9] | |
[10] | |
[11] | |
[12] | |
[13] | |
[14] | |
[15] | |
[16] | |
[17] | |
[18] | |
[19] | |
[20] | |
[21] | |
[22] | |
[23] | |
[24] | |
[25] | |
[26] | |
[27] | |
[28] | |
[29] | |
[30] | |
[31] | |
[32] | |
[33] | |
[34] | |
[35] | |
[36] | |
[37] | |
[38] | |
[39] | |
[40] | |
[41] | |
[42] | |
[43] | |
[44] | |
[45] | |
[46] | |
[47] | |
[48] | |
[49] | |
[50] | |
[51] | |
[52] | |
[53] | |
[54] | |
[55] | |
[56] | |
[57] | |
[58] | |
[59] | |
[60] |